Implications for Primordial Non-Gaussianity {/nl) from weak lensing masses of high-z 

galaxy clusters. 
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The recent weak lensing measurement of the dark matter mass of the high-redshift galaxy cluster 
XMMUJ2235. 3-2557 of (8.5 ± 1.7) x lO" M© at z = 1.4, indicates that, if the cluster is assumed 
to be the result of the collapse of dark matter in a primordial gaussian field in the standard LCDM 
model, then its abundance should be < 2 x 10""^ clusters in the observed area. Here we investigate 
how to boost the probability of XMMUJ2235. 3-2557 in particular resorting to deviations from 
Gaussian initial conditions. We show that this abundance can be boosted by factors > 3 — 10 if the 
non-Gaussianity parameter /^£''' is in the range 150 — 200. This value is comparable to the limit 
for /ivL obtained by current constraints from the CMB. We conclude that mass determination of 
high-redshift, massive clusters can offer a complementary probe of primordial non-gaussianity. 

PACS numbers: cosmology 
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Introduction. — It has been recognized for almost a 
decade that the abundance of the most massive and/or 
high-redshift collapsed objects could be used to con- 
straint the nature of the primordial fluctuation field 
[11) 0, B HL The subject has recently received renewed 
attention E @, S IE 0] possibly sparked by a claimed 
detection of deviations from Gaussianity on CMB maps 
[lo| . Depending on the sign of the non-Gaussian pertur- 
bation, the abundance of rare objects will be enhanced 
or depleted. In [l| we developed the necessary theoreti- 
cal tools to interpret any enhancement (depletion) in the 
abundance of rare-peak objects over the gaussian initial 
conditions case. Working with ratios of non-Gausian over 
the Gaussian case makes the theoretical predictions more 
robust. Later on, Ref. Q generalised the procedure to 
more modern mass-functions and type of non-gaussianity 
including scale dependence. The validity of the analytical 
formulas developed in l| has been recently confirmed by 
detailed N-body numerical simulations with non-gaussian 
initial conditions [7] . These authors have shown that the 
analytical findings in [I| provide an excellent fit to the 
non-Gaussian mass function found in N-body simulations 
with a simple "calibration" procedure. 

Ref. [l3| have recently reported a weak-lensing anal- 
ysis of the z = 1.4 galaxy cluster XMMU J2235.3-2557 
based in HST (ACS) images. Assuming a NFW |ll| 
dark matter profile for the cluster, they estimate a 
projected mass within 1 Mpc of (8.5 ± 1.7) x 10^'' 
M0. Adopting a LCDM cosmology with cosmological 
parameters given by WMAP 5 yr data ([l3l) and 
assuming Gaussian initial conditions they estimate 
that in the surveyed 11 sq. deg. there should be 
0.005 clusters above that mass. Therefore the observed 
cluster is unlikely a the 3(t level. In this Letter we 
explore what level of non-Gaussianity is required to 
boost this abundance by a factor ^ 10 and how this 



relates to the available constraints obtained from the 
CMB. We show that with /^™' in the range 150 - 200 
it is possible significantly enhance the abundance ex- 
pected for such a massive cluster. This value of Jml is 
comparable with current limits from the CMB [13] , [13] ■ 

High Redshift and/or Massive Objects. — While there 
are in principle infinite types of possible deviations from 
Gaussianity, it is common to parameterize these devia- 
tions in terms of the dimensionless parameter /jvl (e.g., 



^ = + /^T'(0'-(0'»- 



(1) 



where $ denotes the primordial Bardeen potential [29| 
and (j) denotes a Gaussian random field. With this con- 
vention a positive value of /^£"' will yield to a positive 
skewness in the density field and an enhancement in the 
number of rare, collapsed objects. 

Although not fully general, this model (called local- 
type) may be considered as the lowest-order terms in 
Taylor expansions of more general fields. Local non- 
gaussianity arises in standard slow roll infiation (al- 
though in this case /^£°' is unnmeasurably small) , and in 



muhi-field models (e.g., [13, 0, [H, |20| ) . For other types 
of non-gaussianity (as we will see below) an "effective" 
Jnl can be defined and related to this model. 

The abundance of rare events (high-redshift and/or 
massive objects) is determined by the form of the high- 
density tail of the primordial density distribution func- 
tion. A probability distribution function (PDF) that pro- 
duces a larger number of > 3f7 peaks than a Gaussian dis- 
tribution will lead to a larger abundance of rare events. 
Since small deviations from Gaussianity have deep im- 
pact on those statistics that probe the tail of the dis- 
tribution (e.g. P, [13]), rare events should be powerful 
probes of primordial non Gaussianity. 
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The non-Gaussianity parameter f^L is effectively a 
"tail enhancement" parameter (c.f., 

As shown in 0, @, 01 when using an analytical approach 
to compute the mass function a robust quantity to use 
is the fractional non-gaussian correction to the Gaussian 
mass function TZ^aiM, z). This quantity was calibrated 
on non-gaussian N-body simulations in [7| . For our pur- 
pose here we want to compute a closely related quantity: 
the non-gaussianity enhancement i.e. ratio of the non- 
gaussian to gaussian abundance of halos above a mass 
threshold Q . As the mass function is exponentially steep 
for rare events here we can safely make the identification 
of the non-gaussianity enhancement with TZng- 

To understand the effect of non-gaussianity on halo 
abundance let us recall that to first order the non- 
gaussianity enhancement is given by 0, 0] : 



7^JVG(M,Z) - 1 + ^3, M 



(2) 



where S^^m denotes the skewness of the density field lin- 
early extrapolated at z = and smoothed on a scale 
R corresponding to the comoving Lagrangian radius of 
the halo of mass M , gm denotes the rms if the -linearly 
extrapolated at z = 0- density field also smoothed on 
the same scale i?; 5'^{zf) = ^5c{zf) and 5c{zf) de- 
notes critical collapse density at the formation redshift 
of the cluster z/. Note that 5c{z) = /:^cD{z = 0)/D{z) 
with D{z) denoting the linear growth factor and Ac is a 
quantity slightly dependent on redshift and on cosmol- 
ogy, which only for an Einstein-de-Sitter Universe is con- 
stant Ac = 1.68. The constant q ~ 0.75 (which we will 
call "barrier factor" ) can be physically understood as the 
effect of non-spherical collapse [2^ lowering the crit- 
ical collapse threshold of a diffusing barrier [26[ see also 
[irt . and has been calibrated on N-body simulations in 
Ref. 't'I. The full expression for TZng is (cf Eqs. 6 and 7 
in Ref. ^]): 



TZng{M,z) = cxp 



3 S3,M 



(3) 



6 ^1 - 5'^S3,m/3 d In <tm 



Let us re-iterate that in principle the enhancement fac- 
tor should be computed by integrating the mass function 
n(Af, z,/7vl) between the minimum and the maximum 
mass and for redshifts above the observed one [4j : 



TZng — 



J n{M, z, fNhjdMdz 
f n{M,zjNL = 0)dMdz 



(4) 



but since the mass function, in the regime we are inter- 
ested in, is exponentially steep, we can identify TZng — 
TZng- 
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FIG. 1: Enhancement factor TZng of the number of rare ob- 
jects for different values of the dark matter mass of the galaxy 
cluster. The lines correspond to different values of /jvl. The 
upper lines are for a collapse redshift of = 2 and the lower 
lines for Zf = 1.4. The shaded area is the range for the 
weak lensing mass estimate of the clusters XMMUJ2235.3- 
2557. Note that for the quoted values of /^£''' it is possible to 
obtain enhancements of order 10 in the cluster number abun- 
dance. This enhancement brings the expected abundance of 
such massive clusters in better agreement with the observa- 
tions. Note that for masses above the estimated central value 
(8.5 X 10^* M0) one expects to find zero such objects in the 
whole sky (one expects 7 objects in the whole sky at the low- 
est value of the mass estimate) which emphasizes the need 
of an enhancement as the one provided by primordial non- 
gaussianity studied here. 



Small deviations from Gaussian initial conditions will 
lead to a non-zero skewness and in particular for local non 
Gaussianity Ss^m = JnT'^Im where denotes the 

skewness produced by /^£''' = 1. Since non-Gaussianity 
comes in the expression for TZng only through the skew- 
ness, the same expression can be used for other types 
of non-Gaussianity such as the equilateral type (see e.g. 
Ref. 0,0] for example of applications). For example, at 



the scales of interest R — 13Mpc/h, S'g'^'^"' 



thus when working on these scales to obtain the same 
non-Gaussian enhancement as a local model, an equilat- 
eral model needs a higher effective /nl- we can make the 
identification /^''^'' = iAf^f^K 

Here we will use the full [l[ expression, corrected for 
the "barrier factor" , for the non-gaussian mass function 
to compute the non-gaussianity enhancement. Note that 
the estimated mass and redshift of XMMUJ2235. 3-2557, 
places it just outside the range where the mass function 
expressions of [l], Q have been directly reliably tested 
with non-Gaussian N-body simulations. Simulations 
seems to indicate that the [l[ expression is a better fit 
than at high masses/redshift and large /atl, this is 
also supported by theoretical considerations Q. 
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Results. — Fig. [T] shows the enhancement factor TZng 
as a function of the mass of the galaxy cluster for different 
values of /^£°' and the redshift of collapse. The shaded 
area shows the error band for the mass determination 
of XMMUJ2235.3-2557 from Ref. [H and the different 
lines have been computed using the [if mass function, 
with the "barrier factor" correction. Ref. 0] show that 
it fits very well the N-body numerical simulations for the 
case of rare peaks, which is the one we are concerned 
with. The solid lines correspond to Jnl — 260, the lower 
one is for a cluster collapse redshift of Zf = 1.4 (i.e. 
assuming that the cluster forms at the observed redshift) 
and the upper one for Zf = 2. The two dashed lines also 
depict the mentioned collapse redshifts but for /nl = 
150. We see that the galaxy cluster abundance can be 
enhanced by a factor up to 10. In the mass range of 
interest, the same enhancement factor can be obtained 
for an equilateral-type non-gaussianity for J^''^*' — 884 
and 510 respectively. 

We should bear in mind that XMMUJ2235. 3-2557 is 
an extremely rare object, sampling the tail of the mass 
function which may not be well known and may be 
strongly affected by cosmology. Using the Q mass func- 
tion we estimate that in the WMAP5 LCDM model (s^l 
one should find 7 galaxy clusters in the whole sky with 
mass greater or equal than the lower mass estimate of 
XMMUJ2235.3-2557 M = 5 x IO^^Mq and z > 1.4 cor- 
responding to a probability of 0.002 for the 11 deg^ of 
the survey. This should be compared with the reported 
number of 0.005 obtained by [13] for a different cosmology 
and different mass function. Thus the effects of cosmol- 
ogy and uncertainty in the mass function can account for 
a factor ~ 2 uncertainty in the predicted halo abundance. 

Note that in all our calculations we have used a conser- 
vative lower limit for the mass of the cluster. If instead 
we use the central or upper value for the mass, using the 
WMAP5 cosmology and the [2^ mass function we ex- 
pect to find zero such clusters in the whole sky, which 
will make our conclusions even stronger. 

The survey area used in Ref. [l^ is 11 deg^, but the 
XMM serendipitus survey in 2006 covered 168 deg^ and 
today covers ~ 400 deg^. Below we report the Ref. [l2| 
numbers and in parenthesis the numbers we obtain. The 
probability of finding XMMUJ2235.3-2557 is thus 0.005 
(0.002) if using 11 deg^; to avoid as much as possible bi- 
ases due to a posteriori statistics one could use 168 deg'^ 
obtaining a probability of 0.07 (0.03), or, as a limiting 
case, even 0.17 (0.07) if using 400 deg^. Note that it is 
likely that there are more clusters as massive in the sur- 
vey area [28| and therefore these numbers are conserva- 
tive. If we use from Fig.[l]the factor 3 tolO enhancement, 
we find that it would help boost the probability to ^ 1 
in the surveyed areas. 

The latest WMAP compilation [13] reports —9 < 

f^' < 111 and -151 < f^^'^ < 253 at 95% confidence, 
reports 27 < /^£"' < 147. The CMB however probes 



much larger scales {R > 120Mpc/h) than those probed 
by clusters such as XMMUJ2235. 3-2557 R - 13Mpc//i: 
a scale-dependent /jvl with k ^ —0.3 can yield an effec- 
tive /jvL on dependence XMMUJ2235. 3-2557 scales that 
is larger than the CMB one by a factor of 3. 

The /^£°' values needed to accomodate the observed 
cluster at z = 1.4 is in the range 150 to 260. This is 
comparable to the limits quoted by Ref. [l^ and [lo| . 

Conclusions. — Accurate masses of high-redshift clus- 
ters are now becoming available through weak lensing 
analysis of deep images. As already discussed in previ- 
ous papers P, Q , their abundance can be used to put 
constraints on primordial non-gaussianity. /j^™' in the 
range 150 — 260 can boost the expected number of mas- 
sive (> 5 X 10^'* M0) high redshift (z > 1.4) clusters by 
factors of 3 to 10. Such large numbers would help make 
clusters like XMMUJ2235. 3-2557 much more probable. 
The scales probed by clusters are smaller than the CMB 
scales, and in principle non-Gaussianity may be scale- 
dependent, making this a complementary approach. 

The adopted error range in the mass determination 
of XMMUJ2235.3-2557 is 100%; even with such a large 
mass uncertainty and considering the pessimistic esti- 
mate of 7 such objects expected in the entire sky with 
a Poisson error of ±2.6, if the entire sky could be cov- 
ered, /^£''' 150 could be detected at > 4(t level. 

RJ and LV acknowledge support of MICINN grant 
AYA2008-03531. LV acknowledges support of FP7- 
PEOPLE-2002IRG4-4-IRG#202182. RJ is supported by 
a FP7-PE0PLE-IRG grant. 



* Email: r aulj , Iverde® astro .princeton.edu' 
[1] Matarrese S., Verde L., Jimenez R., 2000, ApJ, 541, 10 
[2] Robinson J., Gawiser E., Silk J., 2000, ApJ, 532, 1 
[3] Willick J. A., 2000, ApJ, 530, 80 

[4] Verde L., Jimenez R., Kamionkowski M., Matarrese S., 

2001, MNRAS, 325, 412 
[5] Lo Verde M., Miller A., Shandera S., Verde L., 2008, 

JCAP, 4, 14 

[6] Kamionkowski M., Verde L., Jimenez R., 2009, JCAP, 1, 
10 

[7] Grossi M., Verde L., Carbone C, Dolag K., Branchini E., 
lannuzzi F., Matarrese S., Moscardini L., 2009, MNRAS, 
398, 321 

[8] Pillepich, A., Porciani, C, & Hahn, O. 2008, 

arXiv:0811.4176 
[9] Desjacques, V., Seljak, U., & Iliev, I. T. 2009, MNRAS, 

396, 85 

[10] Yadav A. P. S., Wandelt B. D., 2008, PhRvL, 100, 181301 
[11] Navarro J. P., Frenk C. S., White S. D. M., 1997, ApJ, 

490, 493 

[12] Jee M. J., et al., 2009, arXiv, larXiv:0908.3897l 

[13] Komatsu, E., et al. 2009, ApJS, 180, 330 

[14] Coles P., Barrow J. D., 1987, MNRAS, 228, 407 



4 



[15] Verde L., Wang L., Heavens A. F., Kamionkowski, 2000, 

MNRAS, 313, 141 
[16] Komatsu, E., & Spergel, D. N. 2001, PRD, 63, 063002 
[17] Luo X., 1994, ApJ, 427, L71 

[18] Falk T., Rangarajan R., Srednicki M., 1993, ApJ, 403, 
IL 

[19] Gangui A., Lucchin F., Matarrese S., Mollerach S., 1994, 
ApJ, 430, 447 

[20] Fan Z., Bardeen J. M., 1992, University of Washington 

preprint, UW-PT-92-11 
[21] Press W. H., Schechter P., 1974, ApJ, 187, 425 
[22] Fry J. N., 1986, ApJ, 308, L71 
[23] Sheth R., Tormen G.,1999, MNRAS, 308, 119 



[24] Lee J., Shandarin S. F.,1998, ApJ, 500, 14 

[25] Sheth R., Mo H. J., Tormen G., 1 999, a stro-ph/9907024 

[26] Maggiore, M., & Riotto, A. 2009. larXiv: 0903. 1250 

[27] Robertson, B., Kravtsov, A., Tinker, J., & Zentner, A. 

2008, arXiv:0812.3148 
[28] Stanford, S. A., et al. 2006, ApJL, 646, L13 
[29] Which, on sub-Hubble spaces reduces to the usual New- 
tonian potential with a negative sign. 
[30] The cosmological parameter which uncertainty affect 
most the cluster abundance is erg, the value adopted is 
erg = 0.77. 



